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Modeling of Dielectric Barrier Discharge allow for the study of basic plasma processes, and also enable 
analysis and optimization of current technologies using plasmas, as well as predicting the performance of as-yet- 
unbuilt systems for new applications. Modeling of Dielectric Barrier Discharge (DBD) allow for the study of 
basic plasma processes, and also enable analysis and optimization of current technologies using plasmas, as well 
as predicting the performance of as-yet-unbuilt systems for new applications. These discharges are sometimes 
called silent discharges; usually operate at frequencies between 0.05 and 500 kHz. It is known that DBD gener- 
ates non-thermal equilibrium plasma even under atmospheric pressure, which is suitable for ozone production 
[1]. Copious quantities of reactive oxygen species are often generated in atmospheric pressure plasmas. The re- 
active species and amounts generated depend on plasma parameters, such as the rotational, vibrational, excitation 
temperature and electron density [2]. 


This work represents a contribution of modeling of Dielectric Barrier Discharge in oxygen by using 
Comsol Multiphysics software, the powered electrode is driven by sinusoidal high voltage under 50 kHz of fre- 
quency, and the chemical model contains 75 reactions with a various species, where the electron collision with 
the heavy species included in this work are taken from the references [3.4]. The simulation results show an 
ozone molecules density and oxygen density and the others species. 
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A one Dimension symmetric DBD model was developed, 
for the reason of ozone generation, using a sinusoidal high 
voltage. Simulation result show an ozone molecules density and 
oxygen density and the others species. The model contain 75 
reactions with a various species consisting to evaluate the model 
and getting an optimum result. The numerical model used in this 
work was described in abridgement. 


Introduction 


Ozone is a powerful disinfecting and oxidizing agent, and for this reason, 
it is used in a wide range of applications such as treatment of municipal 
and wastewater, food processing, fire restoration, restoration of buildings 
and other objects after floods, etc. When ozone is exposed to organic 
compounds or bacteria, the extra atom of oxygen destroys the 
contaminant by oxidation [1]. 

Ozone generation with Dielectric Barrier Discharge DBD device is an 
interesting search subject for the advantages of DBD’s including the low 
gas flow needed, the homogeneous discharges ignited over comparably 
big area together with different geometries of electrode leading to good 
adaptability [2]. 

Ozone is usually generated between two metallic electrodes, of which 
one electrode is covered with dielectric, using an ac or dc high voltage. 
DBD discharges, which are sometimes called silent discharges, usually 
operate at frequencies between 0.05 and 500 kHz. It is known that DBD 
generates a non-thermal equilibrium plasma even under atmospheric 
pressure, which is suitable for ozone production. [3]. 

Modeling of DBD allow for the study of basic plasma processes, but also 
enable analysis and optimization of current technologies using plasmas, 
as well as predicting the performance of as-yet-unbuilt systems for new 
applications.[5] 

Copious quantities of reactive oxygen species are often generated in 
atmospheric pressure plasmas. The reactive species and amounts 
generated depend on plasma parameters, such as the rotational, 
vibrational, excitation temperature and electron density.[6] 

This proposed model include plasma chemistry, surface reactions, 
density of electrons and heavy species (ions, atom...), electric failed, for 
better understanding DBDs at atmospheric pressure in Oy. 


DESCRIPTION OF THE MODEL 


« Electron Transport 

The continuity equation for the electron density is given by: 

Se tm Ih (1) 
Te = —ueking — V(Dene) (2) 
With: n,: Electron density (1/m*) , R.: Electron rate expression (1/m?.s), 
He: Electron mobility (m?/V .s), E : Electric field (V/s), 
D.: Electron diffusivity (m/s) and e : refers to electron. 
Electron energy is given by: 
OMe + Vp + ETy = Re (3) 
Te = —peEn, — V(Dene) (4) 

With: 
ne: Llectron energy density (V/m*), R-: Energy loss/gain due to inelastic 
collisions (V/(m*.s)), u.: Electron energy mobility (m?/V .s), D.: Electron 
energy diffusivity (m?/s) and e : refers to energy. 

Electron energy, ¢ (V) computed by: 

= = (5) 

The source coefficients in the above equations are determined by the 
stoichiometry of the system and are written using either rate coefficients or 
Townsend coefficients [4]. 

Suppose that there are M reactions which contribute to the growth or decay 
of electron density and P inelastic electron-neutral collisions [4]. 

The electron energy loss is obtained by summing the collisional energy loss 
over all reactions: 
p= en xk) Ny Ne Ag; (6) 

With: P: inelastic clectron-neutral collisions; x;: Mole fraction of the target species 
for reaction j; kj: Rate coefficient for reaction j (m?/s); N,,: Total neutral number 
density (1/m?) and As;: Energy loss from reaction j (V); 


The Townsend coefficients k; depend exponentially on the mean electron 
energy, « - When a Maxwellian EEDF is assumed, the Townsend coefficients 
can be fitted with a function of the form: 

kj = Ashe Fle (7) 


a Efectrostatics Equations 

e Electrostatic field : 

Vege — 10) (8) 
¢ The space charge density : 

p=q(Uh, Zyny — ne) (9) 
e The constitutive relation : 

D = &é,E (10) 

This relation (Y) describe the macroscopic properties of the dielectric coating 
medium (relating the electric displacement D with the electric field E) and the 
applicable material properties. 


* Surface charge accumulation is added to boundaries by way of the 
following boundary condition: 


Di D>) "ps (11) 
where p, is the solution of the following distributed ODE on the boundary: 
ali nj; +nJe (12) 


dt 
Where n./; and n./, are the normal component of the total ion current density 
and the total electron current density respectively on the wall. 


PLASMA CHEMISTRY 


In this part of our work is an aggregation of possible chemical 
reactions in oxygen O, as a gas filled in discharge gap. 


Electron impact reactions : The cross section set for electron 
collisions in O, used in this calculation is shown in Tab.1 [7.8]. 


Reactions and rate coefficient in Tab.2, and the wall reactions in the 


table 3. 
TABLEL Phelps databases 
Reaction Formula Reaction Formula 
1 = O=>e+ 0: 16 | e+Odbis=>e+O2bls 
2 e+ O:=>0+0 17 | e+O2bIs=>e+0+O 
3 e+ O=>e+O2ald 18 | e-O2bls=>2e+ Os" 
4 €+O2ald=>e+ 0: 19 | e+O245=>e+0+0 
5 e+ O:=>e+O2bIs 20 | e+O245=>2e+ O1' 
6 e+O2b1s=>e+ O» 21 
7 e+ O:=>e+0245 22 e+O=>e+Old 
3 e+0245=>e+ O: 23 e+O1d=>e+0 
9 e+ O:=>e+0+0 4 e+O=>e+Ols 
10 | eO=>e0-Old 25 e+O1s=>e+O 
u eO=>e+0-Ols 26 e+O=>2e+0' 
2 e+ O:=>2e+ Or" 27 e+O1d=>e+01s 
13° | e+O2ald=>e+O2ald 28 e+O1d=>2e+O' 
I | eO2ald=>e+0+0 29 e+Ols=>2e+0! 
15 | e+O2ald=>2e+0:' 30 e+ O:=>e+0; 
16 | e-O%bdis=>e+O%bls 3 e+ O=>0+ 0; 
17 | eO2bIsrer0-O 32 ¢ Om>2e10s" 
TABLET. Most important reactions 
Formula Rate constant kt Formula Rate constant kt 
e400) eT haldHO4040 Del 
PO! 001d LBSeLMTe3017 Ona, 404040 Self 2859Te) 
oO L5OMelOTe 65 Oohis-O:+0+0,  231TeD8Te'05 
e000; 3ée-40*Te'0.5 Olbls00+0 — 4.6dde.22*Te'0.5 
0+0;+0+0+ 1e-18*Te'04 OM0;30+0 — Jeld 
Orda +0)+0, Melb emp-$0Te) QO, 400, — 1el6 
OOlldO-e Tel O14#0; +0040 Le.16 
OO Ordre Gell OOO) Set Tay(61Te 
00/3070;  42e7 Old+0 4040 eld 
0:+0+400 De T¥TeBONYOS 15405020; 4Bhel6 
O;+0r400; Tle l6*(Te300)MIS 149, 40,40 4.e18ex(61Te) 
O+0+-700; Te l6*(Te300)05 449 +049 3.3 "ex 30) 
O+0' +00; Abe L6*(Te!300)"05 040,400; Bex 60 Te) 
Or 0X) De TERNS O4040,-40+0; $e Te 083 
OH—Ore Wels O04 Sell Te 
0+0+0;070; — Ge6R(Te300)24 te 40 Lel5 
Oald+ 0,00; Feld erp 200T) ronald +4040 73e-1%er( $90) 
TABLE I. Surface reactions 
Reaction Formula 
£ O:'—0) 
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3 O.'—0; 
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Fig.2. Neutral species densities 


Result Discussion 


When the voltage applied (Fig.2) to the high voltage wire, an electric 
field forms in the gap. Free electrons in the gap will be accelerated 
and may acquire enough energy to cause ionization. The newly 
created electrons rush towards the grounded dielectric circle. An equal 
number of ions are also created since the electrons and ions must be 
created in equal pairs. The ions rush towards the dielectric circle 
where the voltage is being applied with negative polarity. 

From Figure 2, it can be understood that there is a great rise in 
number densities of electron and Oxygen ions the ions and free 
electrons concentration is near the dielectric coatings connected to 
high voltage that means until this time avalanches occur in this area. 
However as the applied voltage increases, a stronger electric field 
forms in the gap between electrodes and any free electron in the gap 
will be accelerated and since the electric field is strong enough they 
acquire enough energy to cause significant ionization activities that 
means avalanches occurs in across the gap uniformly 


Conclusion 


For a large-scale of using ozone in servile industrial 
and medical applications, ozone is mainly prepared by 
non-thermal plasma technologies based on electrical 
discharges. 


This paper develops model to study ozone 
generation by DBD. It simulated with a relatively 
complete physico-chemical model that includes 
elementary plasma processes and chemical reactions 
between discharge species. 


The results of the simulation have shown the time 
evolution of species concentration in gap. 
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